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We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to 
construct an approximation for the exchange-correlation term of the Kohn-Sham approach. The 
resulting exchange-correlation potential is able to capture the features of the strongly-correlated 
regime without breaking the spin or any other symmetry. In particular, it shows "bumps" (or 
barriers) that give rise to charge localization at low densities and that are a well-known key feature 
of the exact Kohn-Sham potential for strongly-correlated systems. Here we illustrate this approach 
for the study of both weakly and strongly correlated model quantum wires, comparing our results 
with those obtained with the configuration interaction method and with the usual Kohn-Sham local 
density approximation. 



I. INTRODUCTION 

In semiconductor nanostructures, the regime of strong 
correlation is reached when the electronic density be- 
comes low enough so that the Coulomb repulsion be- 
comes dominant with respect to the kinetic energy of the 
electrons. From the purely fundamental point of view, 
the study of the strongly-interacting limit in such sys- 
tems is interesting since charge localization, reminiscent 
of the Wigner crystallizatiorti" of the bulk electron gas, is 
expected to occur at low densities. 

A lot of previous theoretical work on Wigner localiza- 
tion in nanostructures has focused on finite-sized quan- 
tum dots (sec, for example, Refs. [2HZ1), and the crossover 
from liquid to localized states in the transport properties 
of the nanostructure has been addressed.'^ ^ In quasi one- 
dimensional nanosystems, signatures of Wigner localiza- 
tion were observed experimentally in one-dimensional 
cleaved-edge overgrowth structures,'^ or in the trans- 
port properties of InSb nanowire quantum-dot systems.l^ 
More recent experimental work clearly identified the for- 
mation of Wigner molecules in a one-dimensional quan- 
tum dot that was capacitively coupled to an atomic 
force microscope probe. Wigner localization has also 
been investigated in other ID systems such as carbon 
nanotubes.l^'l^ (For a review, see Ref. [TB|l . Finally, 
regarding practical applications, Wigner-localized sys- 
tems have been shown to be poten tially useful, e.g., for 
quantum-computing purposes .1^^^^ 

When trying to model electronic strongly-correlated 
systems, however, the commonly employed methodolo- 
gies encounter serious difficulties of different nature. On 
the one hand, the configuration interaction (CI) ap- 
proach, despite being in principle capable of describing 
any correlation regime, is in practice limited to the study 
of small systems with only very few particles due to its 
high computational cost, which scales exponentially with 
the number of particles, N . Such numerical difficulties 
get even worse in the very strongly-correlated limit due 
to the degeneracy of the different quantum states and 



the consequent need of considering larger Hilbert spaces 
in the calculations. Other w ave-function methods like 
Quantum Monte CarlcPEsEl (QMC) and density matrix 
renormalization group (DMRG) ,[2S] which rely to some ex- 
tent on various approximations, can treat systems larger 
than the CI approach, but are still computationally ex- 
pensive and limited to TV < 100. 

The much c heape r Kohn-Sham (KS) density functional 
theory (DFT),'^!!^ which allows to treat thousands of 
electrons, is the method of choice to treat larger quantum 
systems. However, all the currently available approxima- 
tions for the exchange-correlation fu nctio nal fail to de- 
scribe the strongly-correlated regim^^IlSHIZl even at the 
qualitative level. Allowing spin- and spatial-symmetry 
breaking may yield reasonable total energies, without, 
however, capturing the physics of charge localization in 
non-magnetic systems. Moreover, broken symmetry so- 
lutions often yield a wrong characterization of various 
properties and the rigorous KS DFT framework is par- 
tially lost (see, e.g., Refs. 12012312^ . 

KS DFT is, in principle, an exact theory that should be 
able to yield the exact energy and density even in the case 
of strong electronic correlation, without artificially break- 
ing any symmetry. However, when dealing with practical 
KS DFT, one could expect that the non-interacting refer- 
ence system introduced by Kohn and Sham might not be 
the best choice when trying to address systems in which 
the electron-electron interactions play a dominant role. 
For many years, huge efforts have been made in order to 
try to get a better characterization and understanding of 
the properties of the exact Kohn-Sham reference system 
(see e.g. Refs. I20|25|28H44|) . All these works reflected 
the large difficulties encountered when trying to obtain 
adequate approximations to describe strong correlation 
in the exact KS theory.^^ 

An alternative density-functional framework, based 
on the study of the strongly-interacting limit of the 
Hohenberg-Kohn density functional, was presented in 
Ref.Sni In this approach, a reference system with infinite 
correlation between the electrons was considered instead 
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of the non-interacting one of Kohn and Sham. The two 
formalisms can therefore be seen as complementary to 
each other and, indeed, the first results obtained with 
this so-called strictly- correlated- electrons (SCE) DFT, 
presently limited to either ID or spherically-symmetric 
systems, showed its ability to describe systems in the 
extreme strongly-correlated regi me wi th a much better 
accuracy than standard KS DFT.I^SllZI Qn the downside, 
however, SCE DFT requires that one knows a priori that 
the system is in the strong-interaction regime, and it fails 
as soon as the fermionic nature of the electrons plays a 
significant role.^^ Furthermore, the formalism lacks some 
of the appealing properties of the Kohn-Sham approach, 
such as its capability to predict (at least in principle) 
exact ionization energies. Also, crucial concepts widely 
employed in solid state physics and in chemistry, such as 
the Kohn-Sham orbitals and orbital energies, are totally 
absent in SCE DFT. 

Very recently, a new approach that combines the 
advantages of the KS and the SCE DFT formalisms, 
consisting in approximating the Kohn-Sham exchange- 
correlation energy functional with the strong-interaction 
limit of the Hohenberg-Kohn energy density functional, 
has been proposed.!^ Pilot tests of this new "KS SCE" 
framework showed that it is able to capture the features 
of both the weakly and the strongly-correlated regimes 
in semiconductor quantum wires, as well as the so-called 
2kp — >■ Akp crossover occurring in between them, while 
keeping (at least for ID systems) a computational cost 
comparable to the one of standard KS DFT with the 
local-density approximation (LDA). In other words, the 
SCE functional yields a highly non-local approximation 
for the exchange-correlation energy functional, which is 
able to capture key features of strong correlation within 
the KS scheme, without any artificial symmetry break- 
ing. 

The main purpose of this work is to further investigate 
this new KS SCE method, by discussing its exact formal 
properties and, for the prototypical case of (quasi)-lD 
quantum wires, by also performing full CI calculations 
to compare electronic densities, total energies and one- 
electron removal energies in different regimes of correla- 
tion. We find that the KS SCE results are qualitatively 
right at all correlation regimes, representing an impor- 
tant advance for KS DFT. However, while one-electron 
removal energies are quite accurate, total energies and 
ground-state densities are still quantitatively not always 
satisfactory, and therefore we also discuss the construc- 
tion of corrections to KS SCE. In particular, we investi- 
gate here a simple local correction, which, however, turns 
out to give rather disappointing results, suggesting that 
to further improve KS SCE we need semi-local or fully 
non-local density functionals. 

The paper is organized as follows. In the next Sec. |TT] 
we describe the KS SCE approach, illustrating and dis- 
cussing its features beyond what was reported in Ref. |48l 
In Sec.|III|we introduce the quasi-lD systems we have ad- 



the performances of KS SCE with the "exact" CI results, 
with the standard KS LDA method, and discussing KS 
SCE with a simple local correction. Finally, in Sec.[V|we 
draw some conclusions, as well as an outlook for future 
works. 

Hartree (effective) atomic units are used throughout 
the paper. 



II. THEORY AND METHODOLOGY 

A. KS and SCE DFT 

In the formulation of Hohenbcrg and KohrPS the 
ground-state density and energy of a many-electron sys- 
tem are obtained by minimizing the energy density func- 
tional 



E[p]^F[p]+ / drz;e.t(r)p(r) 



(1) 



with respect to the density p{r). In Eq. ([T]) Wcxt(r) is 
the external potential and F[p] is a universal functional 
of the density, defined as the minimum of the internal 
energy (kinetic energy T plus electron-electron repulsion 
Vee) with respect to all the fermionic wave functions ^' 
that yield the density p{r)^ 



F[p] = min(1'|r + Ve 



(2) 



In order to capture the fermionic nature of the elec- 
tronic density, Kohn and SharrPl] introduced the func- 
tional Ts[p] by minimizing the expectation value of T 
alone over all the fermionic wave functions yielding the 
given p{r)^ 



min('J'|r|*), 



(3) 



thus introducing a reference system of non-interacting 
electrons with the same density as the physical, interact- 
ing, one. The remaining part of F[p\, defining the Hartree 
and the exchange-correlation functionals, F[p\ — Ts[p] = 
-EhxcM = E]i[p\ + E^c[p\, is then approximated. The 
minimization of the total energy functional E[p] with re- 
spect to the density yields the well-known single-particle 
Kohn-Sham equation^^H 



1. 



V^-f VKs[p](r) 0.(r) =e,0,(r) , 



(4) 



where WKs(r) = 'yoxt(r) + 5Eyl^c[p\/ 5p{r) = i>cxt[p](r) + 
V h[p\{v) ^ Vxc\p\{^) is the one-body local Kohn-Sham po- 
tential, with vh[p\(j) and Wa;c[p](r) being, respectively, 
the Hartree and the exchange-correlation parts. The so- 
lutions (pi of Eqs. Q are the so-called Kohn-Sham or- 
bitals, which yield the electronic density through the re- 
lation p{y) = |0i(r)p, with the sum running only over 
dressed, and in Sec. IV we present our results, comparing occupied orbitals. Notice that here we work with the 



3 



original, spin-restricted, KS scheme, in which we have 
the same KS potential for spin-up and spin-down elec- 
trons. 

The HK functional of Eq. ^ and the KS functional 
of Eq. ^ can be seen as the particular values at A = 1 
and at A = of a more general functional Fx [p] in which 
the coupling-strength interaction is rescaled with a real 
parameter A, i.e., 



Fx[p] = mmi-^lT + XVe 



(5) 



A well-known exact formula f or the Hartree-exchange- 
correlation functional i^HxcM i^SSM] 



i^HxcW- / {^x[p]\Vee\^x[p])dX= / KeW^^A, (6) 

Jo Jo 

where ^a[p] is the minimizing wave function in Eq. ([s]). 

In the strictly-correlated-electrons DFT (SCE DFT) 
formalism, one considers the strong-interaction limit of 
the Hohenberg-Kohn functional, A — >■ oo, which corre- 
sponds to the functional^^ ■''^ 



yj'.^^M = min(vI.|Ke|*), 



(7) 



i.e., the minimum of the electronic interaction alone over 
all the wave functions yielding the given density p(r). 
This limit has been first studied in the seminal work 
of Seidl and coworker^^^lMI^ and later formalized and 
evaluated exactly in a rigorous mathematical way in 
Refs. i47ii55n57t The functional V^^^[p] also defines a ref- 
erence system complementary to the non-interacting one 
of the Kohn-Sham kinetic energy Tsip], one composed by 
infinitely-correlated electrons, with zero kinetic energy. 
This implies that, analogously as in a set of confined clas- 
sical repulsive charges, which arrange themselves seeking 
for the stable spatial configuration that minimizes their 
interaction energy, in the SCE reference system the po- 
sition of one electron uniquely determines the position of 
the remaining ones, always under the constraint imposed 
by Eq. ([t]) that the density at each point is equal to that 
of the quantum- mechanical system with A = 1, p{r). 

More precisely, the functional V^^^ [p] is constructecP^ 
by considering that the admissible configurations 
of N electrons in d dimensions are restricted to 
a d— dimensional subspace flo of the full classical 
iVd— dimensional configuration space. A generic point 
of fio has the form 



Roa(s) = (fl(s), ,fjv(s)). 



(8) 



where s is a d-dimensional vector that determines the 
position of, say, electron "1", and fi(s) {i — 1,...,N), 
with fi(s) — s, are the so-called co-motion functions, 
which determine the position of the «-th electron as a 
function of s. The co-motion functions ar e implicit n on- 
local functionals of the given density p{r) HSMElM] ^^^^^ 
solution of a set of differential equations that ensure the 



invariance of p under the coordinate transformation s — >■ 
fj(s), i.e., 

p(f,(s))df,(s) =p(s)ds, (9) 

or, equivalently, that the probability of finding the elec- 
tron i at fi{s) is equal to that of finding the electron 
"1" at s. At the same time, the ii{s) must satisfy group 
propertie s tha t ensure the indistinguishability of the N 
electrons .'^^'^ 

The functional Vj,^^ [p] can then be written in terms 
of the co-motion functions a d^^ * ^^ * 



N 



(10) 



just as Ts[p] is written in terms of the Kohn-Sham or- 
bitals ^i(r). The equivalence of the two expressions for 
V^^^^Ip] in Eq. (fTol) has been proven in Ref. 59', 



Since in the SCE system the position of one electron 
determines all the other — 1 relative positions, the 
net repulsion felt by an electron at position r due to 
the other A^ — 1 electrons becomes a function of r it- 
self. For a given density po(r), this effect can be ex- 
actly transforme d^^ ' ' ^^ into a local one-body effective 
external potential t'sCE[po](r) that compensates the to- 
tal Coulomb force on each electron when all the particles 
are at their respective positions fJpoKr), i.e., such thalP^ 



V-eWW^g^fffi., (11) 

In terms of the classical-charge analogue, wscE[po](r) can 
thus be seen as an external potential for which the total 
classical potential energy 

N-l N ^ N 

-Epot(i-i,...rjv) = ^ , ^ ^ I +X!"scE[po](r») 

i=l j=i+l I * J I i=l 

(12) 

is minimum when the electronic positions reside on the 
subset R^o, i.e., when r,; = fi[po](r) or, equivalently, 
when the associated density at each point is equal to 
Po(r). For an arbitrary density p(r), the potential-energy 
density functional defined as 

iSpSCE ^ ^^SCE ^ I ^^^^ (r)p(r)dr (13) 

will then satisfy the stationarity condition 
6El^,^[p]/Spir)\^^^^ = 0, i.e., we wiU have that 



Sp{r) 



= -«sce[po](i") 



(14) 



P=Po 
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Notice that Eq. ( 14 1 involves the functional derivative of 



a highly non-local implicit functional of the density, de- 
fined by Eqs. (|9|-([To|. This, however turns out to reduce 



to a local one-body potential that can be easily calculated 
from the integration of Eq. ( 11 ) once the co- motion func- 
tions are obtained via Eq. (|9|. This shortcut to compute 
the functional derivative of V^^^ [p] is extremely powerful 
for including strong-correlation in the KS formalism .I^Sl 



B. Zeroth-order KS-SCE approach 



Equations ( 11 1 and ( 14 ) show how the effects of strong 
correlation, captured by the limit A — > oo of F\[p] and 
rigorously represented by the highly non-local functional 
Vgg^^ [p] , are exactly transferred into the one-body poten- 
tial wsceM- The KS SCE approach to zeroth order con- 
sists in using this property to approximate the Hartree- 
exchange-correlation term of the Kohn-Sham potential 



SEukcIp] 
Sp{r) 



vscE [p] (r) , -DscE [p] (r) = - wsce [p] (r) 



(15) 

Notice that we have defined vsce[p]{i^) = —vsce[p]{t^) , as 
here we seek an effective potential for KS theory, which 
corresponds to the net electron-electron repulsion acting 
on an electron at position r, while the effective poten- 
tial for the SCE system of Eq. (11) compensates the net 
electron-electron repulsion. 

More rigorously, by considering the A — » oo expansion 
of the integrand of Eq. Q one obtain^sHSSlSSl 



A— >c 



SCEr 



o{x-n, 



(16) 



where the acronym "ZPE" stands for "zero-point en- 
ergy", and p >5/4 - see Ref. 158 for further details. By 
inserting the expansion of Eq. ( 16 1 into Eq. (|6| one ob- 



tains an approximation for i?Hxc [p\ , 



(17) 



We consider here only the first term, corresponding to a 
zeroth-order expansion around A = oo, i.e., -Ehxc[p] ~ 
V^^^Ip], which yields Eq. (151 for the corresponding 
functional derivatives. 

Taking into account the definition of the functional 
Vee^^[p], Eq. ([t]), the zeroth-order KS SCE is equivalent 
to approximate the minimization over '5 in the HK func- 
tional of Eq. ([2| as 



min(*|T-|-t/ee|*) 



min(«'|r|1') 



min(*|Vee|*) 



(18) 



The KS SCE approach thus treats both the kinetic energy 
and the electron-electron repulsion on the same footing, 
combining the advantages of both KS and SCE DFT and 



therefore allowing one to address both the weakly- and 
the strongly-interacting regime, as well as the crossover 
between them.l^ Indeed, from the scaling propertied^ of 
the functionals F[p], Ts[p] and V^^^[p] it derives that the 



approximation of Eq. ( 18 1 becomes accurate both in the 



weak- and in the strong-interaction limits, while probably 
less precise in between. To use the scaling relations'^ one 
defines, for electrons in D dimensions, a scaled density 



P7W = 7^P(7I") 



7 > 0. 



We then havePlSl 



F[p^] = j^F,/^[p], 



(19) 
(20) 
(21) 



where meand^Sl that the Coulomb coupling con- 

stant A in Fx[p] of Eq. ([5| has been set equal to I/7. We 
then see that both sides of Eq. (18) tend to Ts[p.y] when 



7 — > 00 (high-density or weak-interaction limit) and to 
Vgg'^^[p^] when 7 — (low-density or strong-interaction 
limit). 

Standard KS DFT emphasizes the non-interacting 
shell structure, properly described through the functional 
Ts[p], but it misses the features of strong correlation. 
SCE DFT, on the contrary, is biased towards localized 
"Wigner-like" structures in the density, accurately de- 
scribed by V^^^[p], missing the fermionic shell struc- 
ture. Many interesting systems lie in between the weakly 
and the strongly interacting limits, and their complex 
behavior arises precisely from the competition between 
the fermionic structure embodied in the kinetic energy 
and correlation effects due to the electron-electron repul- 
sion. By implementing the exact 'DscE[p](r) potential in 
the Kohn-Sham scheme, we thus let these two factors 
compete in a self-consistent procedure.!^ 

One should also notice that while the KS SCE ap- 
proach does not use explicitly the Hartree functional, the 
correct electrostatics is still captured, since Vgg^^[/3] is 
the classical electrostatic minimum in the given density 
p. Moreover, the potential uscE[p](r) stems from a wave 
function (the SCE on^^^^ and is therefore completely 
self-interaction free. 

Finally, another neat property of the zeroth-order KS 
SCE approach is that it always yields a lower bound to 
the exact ground-state energy Eg = E[po] , where po is the 
exact ground-state density. In fact, for any given p the 
right-hand side of Eq. ( 18 ) is always less or equal than 



the left-hand side, as the minimum of a sum is always 
larger than the sum of the minima. As a consequence, 
for p = pq we have the inequality 

E[po] = F[po]+ J pov^^t > T,[po]+K!°^[Po]+ J Po«cxt , 

(22) 

which becomes even stronger when ones minimizes the 
functional on the right-hand-side with respect to the den- 
sity within the self-consistent zeroth-order KS SCE pro- 
cedure. It should be noted that this property implies an 
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important difference with respect to the variational wave- 
function methods (such as HF, CI, QMC and DMRG), 
which, instead, provide an upper bound to the exact 
ground-state energy. 



C. Local correction to zeroth-order KS SCE 

As preliminary found in Ref. |3H! and further shown 
in Sec. |IV[ the zeroth-order KS SCE yields results that 
are qualitatively correct in the strong-correlation regime 
(representing a significative conceptual advance for KS 
DFT), but still with quantitative errors, which become 
smaller and smaller as correlation increases. An impor- 



tant issue is thus to add corrections to Eq. (18). One 
can, more generally, decompose F[p\ as 

F[p] - Ts [p] + ^ [p] + [p] + [P] , (23) 
where Tc[p] (kinetic correlation energy) is 

TM^{m\T\m)-Ts[p], (24) 

i.e., the difference between the true kinetic energy and 
the Kohn-Sham one, and V^J^p] (electron-electron decor- 
relation energy) is 



KeW-(*[p]|Ke|*W)-Kf''[p] 



(25) 



i.e., the difference between the true expectation of Vee 
and the SCE value. A "first-order" approximation for 
Tc[p] + V/Jp] can be obtained from Eq. ([l7]), 

rc[p]+KlW«2Kf^[p], (26) 

and can be, in principle, included exactly using the for- 
malism developed in Ref. i58j but other approximations, 
e.g. in the spirit of Ref. [6T1 can also be constructed. 

Here we consider an even simpler approximation, 
Tc[p\ + V/e[p] w £'lcM, where £^lc[p] is a local term 
that includes, at each point of space r, the corresponding 
correction for a uniform electron gas with the same local 
density p(r), i.e., 



Ehc [p] 



pir)[tMr)) + viipir))]dr. (27) 



In Eq. (27) tc{p) and vf^{p) are the kinetic correlation 
energy and the electron-electron decorrelation energy per 
particle of an electron gas with uniform density p, corre- 
sponding to 

Up) + «fe(p) = e^c(p) - escE(p), (28) 

where exc{p) is the usual electron-gas exchange- 
correlation energy and escE(p) is the indirect part (ex- 
pectation of Vee minus the Hartree energy) of the SCE 
interaction energy per electron of the uniform electron 
gas with density p. This correction makes the approxi- 
mate internal energy functional 

SlcM (29) 



F[p\=Ts[p] 



become exact in the limit of uniform density, similarly to 
what the LDA functional does in standard KS DFT. 



III. MODEL AND DETAILS OF THE 
CALCULATIONS 

We consider N electrons in the quasi-one-dimensional 
(QID) model quantum wire of Refs. 1271621 



^ iv ^9 iv — J. iv A'^ 

m 

in which the effective electron-electron interaction is ob- 
tained by integrating the Coulomb repulsion on the lat- 
eral degrees of freedom,'^^] and is given by 



The parameter h fixes the thickness of the wire, set to 
6 — 0.1 throughout this study, and erfc(x) is the comple- 
mentary error function. The interaction Wb{x) has a long- 
range coulombic tail, Wbix — > oo) = 1/x, and is finite at 
the origin, where it has a cusp. As in Ref. 1271 we con- 
sider an external harmonic confinement Wcxt(a;) = ^uP'x'^ 
in the direction of motion of the electrons. The wire 
can be characterized by an effective confinement-length 
parameter L such that 
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Z2' 



/ \ 1-2 2 
Vcyit[x) ^ -LU X . 



Zeroth-order KS SCE 



The co-motion functions fi{x) can be constructed by 
integrating Eqs. ^ for a given density p(a;),'^2lllIII! choos- 
ing boundary conditions that make the density between 
two adjacent strictly-correlated positions always inte- 
grate to 1 (total suppression of fluctuations),'^ 



p{x) dx = 1, 



(32) 



and ensuring that the fi{x) satisfy the required group 
properties. ^^''^^'^ This yields 



/i(^) ^ ^ AT-i 



Ne^[Ne{x) + I ^ I] 



X < ajv+i-i 



Ne^[Neix)+i-l- N] X>aN+l^^, 

(33) 

where the function Ne{x) is defined as 

Ne{x) = / p{x')dx', (34) 



and flfc = ^{k). Equation (11) becomes in this case 

N 

iscE[p]{x) ^^Wb{\x - fi{x)\)sgn{x - Mx)). (35) 

We then solve self-consistently the Kohn-Sham equations 
^ with the KS potential VKsix) = Vcxt{x) + vsce[p]{x), 
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where vsce [p] (x) is obtained by integrating Eq. ( 35 1 with 
the boundary condition wsCEMda;! — >■ oo) = 0. As said, 
we work in the spin-restricted KS framework, in which 
each spatial orbital is doubly occupied. 



B. The configuration interaction method (CI) 

In the configuration interaction calculations, the full 
many-body wavefunction is expanded as a linear combi- 
nation of Slater determinants, constructed with the non- 
interacting harmonic oscillator orbitals. A matrix repre- 
sentation of the Hamiltonian in this basis is then numer- 
ically diagonalized to find the eigenstates of the system. 
The number of possible ways to place N particles in a 
given set of orbitals increases rapidly as a function of 
such that only small particle numbers are tractable. 
Also, the stronger the interaction, the more basis or- 
bitals are generally required to obtain a good approxi- 
mation. For the present physical system, about 20-40 
orbitals were needed to get converged solutions, which 
resulted in Hilbert space dimensions in the range 10^- 
10®. For a more detailed description of the method, see 
e.g. Refs. [55IM1 



C. KS LDA 

We have performed Kohn-Sham LDA calculations us- 
ing the exchange-correlation energy per particle Cxc — 
Ex + Ec for a ID homogeneous electron gas with the 
renormalized Coulomb interaction Wbix), as detailed in 
Ref. \27i The exchange term ex is given by 



+ C30 



dq 
2^ 



M<i) [Saiq) - 1] 



(36) 



where Vh{q) is the Fourier transform of the interaction po- 
tential, 5*0 (9) is the non-interacting static structure fac- 
tor, and Ts = ^^To increase the numerical stability, 
we have interpolated between the Taylor expansions of 
^xi'f's) at small and large up to order 14. For the cor- 
relation term we have used the results of Casula et al.^ 
who have parametrized their QMC data as 



A + Br]' + Cr^^ 



lii{l + Drs+Er]^), (37) 



where the different parameters are given in Table IV of 
Ref. ESI for several values of b. 
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FIG. 1: (color online) The indirect SCE energy escE(?'s) for 
the ID gas [interaction of Eg. (|31| l and fe = 0.1] is compared 
to the parametrized QMC datcP^for the exchange-correlation 
energy txdrs)- For the SCE energy we show both our numer- 
ical results and the fitting function of Eqs. (38l-(39l. 



have then evaluated the limit A'^ — cxi at fixed density p 
to obtain the bulk value. The details of this calculation 
are reported in Appendix [A] 

In Fig. [1] we show our numerical results for 6 = 0.1 
compared to the parametrizecPH QMC results for the 



exchange-correlation energy e2;c(rs) of Eqs. (36)-(37|. We 
see that, as it should be, escE(''s) < ^xcifs) everywhere. 
For large we find that the SCE data are very close 
to the QMC parametrization, with differences of the or- 
der of ^ 0.1%. Notice also that at = we have 

escE(O) = exc{0) — ^xiO) — This is due to the fact 

that in the — limit at fixed b the first-order per- 
turbation to the non-interacting gas is just a constant, 
so that every normalized wave-function yields the same 
result for the leading term. We have parametrized our 
data for esCE(^s) as 



escE(p) ^ pq{2bp), 



with 



q{x) = Ai In 



1 + azx + 02^2 



(38) 



(39) 



and Ai = 
2.025166778, 



ai 



0.9924534, aa 

a2v^ 
2Ai ■ 



0.3 



= 1.55176743, 03 = 
This fit is valid for all 



values of 6, since the scaling of Eq. ( 38 ) is exact for the 
SCE energy. The fitting function is also shown in Fig. [T] 
for the case b — 0.1. 



D. KS SCE with local correction 



IV. RESULTS 



We have obtained the indirect SCE energy per electron 
escE(p) needed in Eq. ( 28 ) by first computing the indirect 
egQg(p, A^) for a ID droplet with A'' electrons, uniform 
density p and radius R — |^ , as described in Ref. |56l We 



Figure [2] shows the electron densities for A^ = 4 and 
different effective confinement lengths L = 2lu^^^'^ ob- 
tained with the KS SCE, the CI and the KS LDA ap- 
proaches. One can see that the three methods show 
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qualitative agreement in the weakly-correlated regime, 
represented here in panel (a) by the case L = 1. The 
densities have N/2 peaks, given by the Friedel-like oscil- 
lations with wave number 2k'^p ^ where fc^ ~ T^p/S is the 
effective Fermi wavenumber, determined by the average 
density in the bulk of the trap p. 

As the confinement length of the wire increases, the 
interactions start to become dominant and, whereas the 
KS SCE and the CI results are still in qualitative agree- 
ment, the LDA clearly provides a physical wrong descrip- 
tion of the system. Indeed, one can see from panel (b) 
that whereas the densities obtained from the KS SCE 
and the CI methods develop a four-peak structure, cor- 
responding to charge localization and indicating that the 
system enters the crossover between the weakly and the 
strongly correlated regimes (the 2kp — >■ Akp crossover), 
the KS LDA yields a flat density. This is a typical er- 
ror of local and semilocal density functionals that shows 
up also in bond breaking (yielding wrong molecular dis- 
sociation curves) and in systems close to the Mott in- 
sulating regime. In such cases, better total energies are 
obtained by using spin-dependent functionals and allow- 
ing symmetry breaking. This, however, does not yield a 
satisfactory physical description of such systems, missing 
many key features and giving a wrong characterization 
of several properties (see, e.g., Refs. I20l23l26p . 

When the system becomes even more strongly- 
correlated, here represented by L = 70, the KS SCE gets 
closer to the CI result, with densities that clearly present 
N peaks, corresponding to charge localization. The KS 
LDA density is now very delocalized and almost flat in 
the scale of Fig. [2] In order to obtain charge localization 
within the restricted KS scheme, the self-consistent KS 
potential must build "bumps" (or barriers) between the 
electrons. These barriers are a very non-local effect and 
are known to be a key property of the exact Kohn-Sham 
potential, as discussed in Refs. 1^ and 1^ 

In Fig. [3] we show that the self-consistent KS SCE 
scheme builds, indeed, the above-mentioned barriers in 
the corresponding Kohn-Sham potentials, which we plot 
together with the corresponding densities for = 4 and 
AT = 5 for L = 70. One can see that each of the A^ 
peaks in the density corresponds to a minimum in the 
KS potential, which is separated from the neighboring 
ones by barriers or "bumps", at whose maxima the KS 
potential has a discontinuous (but finite'^ first deriva- 
tive. The number of such barriers is thus equal to A^ — 1, 
and they become more pronounced with increasing corre- 
lation, enhancing the corresponding charge localization. 
Notice that the discontinuous first derivative of the KS 
SCE potential at the barrier maxima is a feature due to 
the classical nature of the SCE potential, and it is not 
expected to appear in the exact KS potential (indeed, it 
does not appear in any of the available calculations of 
the "exact" KS potential obtained by inversion). 

It is also interesting to make a connection between 
our results and the recent work on the KS exchange- 
correlation potential for the ID Hubbard chains. ^4=^^^*^ 
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FIG. 2: (color online) Electron densities for N = A and 
L = 1, 15 and 70, obtained with the KS SCE, CI and LDA 
approaches. The results are given in units of the effective 
confinement length L = 2u)' 



^1/2 



In particular, VieirgP^ has shown that the exact 
exchange-correlation potential for a ID Hubbard chain 
with hopping parameter t and on-site interaction U , ob- 
tained by inversion from the exact many-body solution, 
always oscillates with frequency Akp, while the density 
oscillations undergo a 2kp — Akp crossover with in- 
creasing U jt. The crossover in the density is thus due 
to the increase in the amplitude of the oscillations of the 
xc potential. In Fig. |4] we show the KS SCE exchange- 
correlation potentials for A^ = 4 electrons in the weakly 
(L — 2) and strongly (L = 70) correlated regimes. We see 
that the KS SCE self-consistent results are in qualitative 
agreement with the findings of Vieira.'^the oscillations in 
the xc potential have essentially a frequency Akp also in 
the weakly-correlated case, with amplitude that increases 



with increasing L [due to the scaling of Eqs. ( 19 )-( 21 1 the 



parameter L plays here a role similar to \J jt for the Hub- 
bard chain] . In the two lower panels of the same figure 
we also further clarify the 2kp — >■ ^kp crossover in the 
KS framework: we see that the Akp regime in the density 
oscillations occurs when the barriers in the total KS po- 
tential (due to the large oscillations of the xc potential) 
are large enough to create classically-forbidden regions 
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FIG. 3: (color online) Self-consistent Kohn-Sham potentials 
obtained with the KS-SCE method for iV = 4 and TV = 5, 
with effective confinement length L = 70 (blue solid lines). 
The corresponding densities are also shown (red dotted lines). 
Notice that for the sake of clarity only the results for x > 
are shown. The results are given in units of the effective 
confinement length L = 2uj~^^^ . 



inside the trap for the occupied KS orbitals. 

In Table |T] we report the total energies obtained with 
the three approaches, KS SCE, CI and KS LDA, for dif- 
ferent values of the parameters L and N. It can be seen 
that in the weakly-correlated regime, represented here by 
L = 1 and 2, the error made by the KS SCE approach is 
larger than the one corresponding to the KS LDA. The 
results also clearly show that, as previously discussed, 
KS SCE is always a lower bound to the total energy. As 
the system becomes more correlated, the results obtained 
with the KS SCE and the CI approaches become closer 
to each other, whereas the value given by the KS LDA is 
less accurate, as one could have inferred from the corre- 
sponding densities shown in panels b) and c) of Fig. [2] 

In the exact Kohn-Sham theory, the highest occu- 
pied KS eigenvalue is equal to minus the ex act ch em- 
ical potential from the electron-deficient sidej ^^ * ^" -^ i.e., 
fx~ = En^i — Ej^. In Table [lT| we compare the highest oc- 
cupied KS eigenvalue obtained with the KS SCE and the 
KS LDA approaches with the values of En — Em-i cal- 
culated from the total energies given by the CI method, 
corresponding to the same values of N and L given in 
Table |l] One can see that in this case the KS SCE gives 



FIG. 4: (color online) Top panel: the self-consistent KS SCE 
exchange-correlation (xc) potential for A'' = 4 at weak corre- 
lation (L = 2) and strong correlation (L = 70). In the inset, 
the oscillating part of the xc potential at L = 2 is zoomed 
in. Middle panel: the total self-consistent KS SCE potential 
(blue, solid line), the corresponding density (red dotted line), 
and the two occupied KS eigenvalues (green dashed horizontal 
lines) for the weakly-correlated L = 2 wire. In this case, we 
see that in the KS system there are no classically-forbidden 
regions inside the trap. Bottom panel: the same as in the 
middle panel for the strongly-correlated L = 70 wire. In this 
case, we clearly see the classically-forbidden regions inside the 
trap created by the barriers in the KS SCE potential. The 
results are given in arbitrary units. 



N L KS SCE CI KS LDA 



2 


2 


1.81 


2.49 


2.59 


2 


15 


0.0942 


0.106 


0.130 


2 


70 


0.0112 


0.0115 


0.0182 


4 


1 


25.08 


28.42 


28.57 


4 


2 


8.46 


10.60 


10.68 


4 


15 


0.491 


0.541 


0.580 


4 


70 


0.0602 


0.0629 


0.0771 


5 


15 


0.787 


0.871 


0.915 


5 


70 


0.099 


0.102 


0.121 



TABLE I: Comparison of the total energies obtained with 
the KS SCE, CI and KS LDA approaches for different values 
of the particle number A'^ and effective-confinement length 
i = 2a;-i/2. 
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N 


L 


KS SCE 


CI 


KS LDA 


2 


2 


1.65 


1.99 


2.56 


2 


15 


0.104 


0.097 


0.263 


2 


70 


0.0126 


0.0111 


0.04087 


4 


1 


11.26 


11.86 


12.56 


4 


2 


4.08 


4.65 


5.02 


4 


15 


0.248 


0.256 


0.453 


4 


70 


0.0318 


0.0304 


0.06909 


5 


15 


0.325 


0.330 


0.539 


5 


70 


0.0408 


0.0391 


0.08172 



TABLE II: For the same systems of Table |I] we compare the 
highest occupied KS eigenvalues obtained from KS SCE and 
KS LDA with the full CI values of En - En-i- 



good results also in the weakly-correlated regime. In 
the strongly-correlated limit, the KS SCE and the CI re- 
sults show an agreement similar to that observed in the 
corresponding total energies. KS LDA, as usual, yields 
too high eigenvalues, due to the too fast decay of the 
exchange-correlation potential for \x\ — > oo. 

As mentioned earlier, the numerical cost of the CI 
method increases exponentially with the number of par- 
ticles, and this limitation becomes stronger as the cor- 
relations become dominant. In the calculations reported 
above, for the 5-electron case with L =70 we diagonal- 
ized a matrix where the eigenvectors had a dimension of 
about 3.5 X 10^. While it is technically possible to treat 
larger matrices, the rapid growth of the basis size still 
efficiently limits the number of particles one can han- 
dle. (For iV = 6 electrons, using the same basis orbitals, 
the corresponding dimension is roughly 2.6 x 10^.) The 
KS SCE method, on the contrary, has a numerical cost 
(in ID) comparable to the one of KS LDA, therefore al- 
lowing to study strongly-correlated systems with much 
larger particle numbers. In Fig. [5] we show the electron 
densities and corresponding KS potentials obtained with 
the KS SCE method for TV = 8, 16, and 32, for different 
values of L: in panels a) and b) we see how, at fixed 
number of particles = 8, the bumps in the KS poten- 
tial and the amplitude of the density oscillations become 
larger with increasing L. For fixed effective confinement 
length L = 150, we see from panels b), c) and d) how in- 
creasing the particle number N leads to less pronounced 
features of strong correlation, according to the scaling of 
Eqs. (pi- 1211). 



Finally, we have tested the local correction to the 
zeroth-order KS SCE discussed in Sees. Ill Cl and HITdI as 
we see in the case N = 2 and L = 20 reported in Fig. |6j 
the results for the self-consistent densities are very disap- 
pointing, laying in between the KS SCE and the standard 
KS LDA values. This is due to the fact that, similarly 
to the standard KS LDA case, this simple local correc- 
tion cannot capture the physics of the intermediate and 
strong-correlation regime, so that its inclusion worsens 
the results of KS SCE. In future work we will explore 
semi-local and fully non-local corrections to KS SCE. 




C/2 
> 



X 



FIG. 5: (color online) Electron density and corresponding KS 
SCE potential for different particle numbers A'^ and effective 
confinement lengths L. As in Fig.jS] only the results for a; > 
are shown. The results are given in arbitrary units. 




FIG. 6: (color online) Electron density for the case N = 2 
and L = 20. The "exact" CI result is compared with the KS 
LDA the KS SCE a nd the KS SCE with local correction of 
Sees. [ncl and [IiTd] (KS SCE-I-LDA) results. The results are 
given in units of the effective confinement length L = 2lu~^^'^. 
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V. CONCLUSIONS AND PERSPECTIVES 

We have used the exact strong-interaction hniit of 
the Hohenberg-Kohn functional to approximate the 
exchange-correlation energy and potential of Kohn-Sham 
DFT. By means of this so-called KS SCE approach, we 
have addressed quasi-one-dimensional quantum wires in 
the weak, intermediate and strong regime of correlations, 
comparing the results with those obtained by using the 
configuration interaction and the KS local density ap- 
proximation. In the weakly-correlated regime, the three 
approaches give qualitatively similar results, with elec- 
tronic densities showing N jl peaks, associated with the 
double occupancy of the single-particle levels that domi- 
nate the system. In this regime, KS LDA performs over- 
all better than KS SCE. As correlations become domi- 
nant, the KS SCE and the CI densities start to develop 
additional maxima, corresponding to charge-density lo- 
calization, whereas the KS LDA provides a qualitatively 
wrong description of the system, yielding a very flat, delo- 
calized, density. We have also investigated a simple local 
correction to KS SCE, which, however, gives very dis- 
appointing results. In future works we will thus explore 
semi-local and fully non-local corrections to KS SCE. 

The Kohn-Sham potential of the KS SCE approach 
shows "bumps" that are responsible for the charge local- 
ization and are a well-known feature of the exact Kohn- 
Sham potential of strongly-correlated systems. More- 
over, the associated KS SCE exchange-correlation poten- 
tial shows the right asymptotic behaviour, since it is self- 
interaction fr ee as it is constructed from a wave function 
(the SCE onePsMIIil). This way, KS SCE is able to also 
give rather accurate chemical potentials. Notice that, as 
shown by studies of one-dimensional Hubbard chains, the 
2fc^ Akp crossover in the density is a ve ry cha llenging 
task for KS DFT for non-magnetic systems.E^^ The fact 
that KS SCE is able to capture this crossover is thus a 
very remarkable and promising feature. 

Crucial for future applications is calculating Vgg^^[p] 
and wscE[p](r) also for general two- and three- 
dimensional systems. An enticing route towards this goal 
involves the mass-transportation-theory reformulation of 
the SCE functional,'^ in which Vgg^^[(o] is given by the 
maximum of the Kantorovich dual problem, 

/ u{r)p{r)dr : ^u{r^) < ^ XI uTrT 

i=l 1=1 j>i ' * 

where u{r) = vsce[p]{''c) + C , and C is a constant This 
is a maximization under linear constraints that yields in 
one shot the functional and its functional derivative. Al- 
though the number of linear constraints is infinite, this 
formulation may lead to approximate but accurate ap- 
proaches to the construction of Vgg*-^^[p] and wscE[p](r), 
as very recently shown by Mendl and Lin.*^ 
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Appendix A: SCE for the uniform QID electron gas 

Following Ref. [56l we have computed the SCE indirect 
Coulomb interaction energy per electron ^sc^iP'-^) '-'^ 
a ID droplet of uniform density p and radius R — 
where N is the number of electrons, 



N 



2bp 



where 



^SCE 



N 

(x, ^) = ^ - z)e''/^'erfc 



(Al) 



(A2) 



is the rescaled SCE energy of the droplePs! and the sec- 
ond term in the right- hand-side of Eq. ( Al ) is its Hartree 
energy, with 



ui(x) 



sin k 



^ Ei{k^x^)dk, (A3) 



and 



Ei{x) = 



t 



-dt. 



Since the function ui{x) is numerically unstable, we have 
interpolated between its small- a; expansion through or- 
ders O(x^), 



4 7-1 
TTX TTX 1 



1 



- ^-^ + ^7r'^'a:-7rlog(x)-^7rV'^"^ y^j, 



with (I 
through orders 0{x~^^), 



0.036489974, and its large- a; expansion 



7r3/2 647r tt^/^ 

1209600x15 " 14189175x14 ^ 131040xi3 ~ 
167r 7r3/2 167r n^^^ 2n 

405405x12 ^ 15840x11 ~ 51975xio ^ 2160x9 ~ 945x8 

7r3/2 47r Tr^_^ tt^ ^ tt^ 
"^336x7 315x6 ^ 60x5 15a? ^ 12^ 3^ 

switching between them at x = 0.584756. 

We have then evaluated numerically the limit N ^ oo 
of Eq. (Al) at fixed p. We have found that the conver- 
gence is reasonably fast: for example, taking N = 10^ 
yields results with a relative accuracy of 10^^ Our nu- 
merical results have been fitted with the function q(x) of 

Eqs. (Issl-pgl). 
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